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ABSTRACT 

The Module for Automatic Dock and Detumble (MADD) is an automated 
device for bringing a passive, tumbling space base under control in an 
orbital rescue, situation. The conceptual design of -such a device 
resulted from a consideration of tumbling motion analyses, and mission- 
constraints. Specific topics of investigation include orbit and. 
attitude dynamics and detumble profiles. Position and attitude control 
systems for the various phases of operation were developed. Dynamic 
motion of a passive vehicle with MADD attached is considered as an 
^example application -and to determine control requirements. Since. time 
is a critical factor in rescue operations, it is essential to execute 
the detumbling maneuver in a minimum of time. Optimization of the 
MADD thrusting sequence has also been investigated. Results indicate 
the control torque must be directed opposite to the angular momentum 
vector for the assumptions used here. 
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Nomenclature 


principal moments of inertia of MSS 

3x3 direction cosine matrix defining the actual 
"orientation of MADD with respect to the command 
orientation (error matrix) 

.^angular velocities in a right hand system 

3 x 3 direction cosine matrix defining the. attitude 
of MADD with respect to the moving coordinate f-rame 

3x3 direction cosine matrix defining the- command 
attitude for MADD with respect to the moving- coordinate 
frame. 

unit orthogonal vectors attached to the body fixed 
axes of MADD defining its coordinate frame, 

-(i = 1,2,3) 


components of the eigenvector b in the body roll , 
pitch, and yaw axes, respectively, (i = x,y,z) 

unit orthogonal vectors defining the command attitude 
of MADD, (i = 1,2,3) 

moments of inertia of MADD, (i = x,y,z) 

-gain terms associated with the control law about 
the body roll, pitch, and yaw axes K = [I. L(max)/ 

I (max) ] (i = x,y , z) A± ' 1 

rate gain term associated with the control law, 

1/w (max) 
a 

values of control torque, (i = x,y,z) 

unit orthogonal vectors defining the moving (reference ) 
coordinate frame, (i = 1,2,3) 

angle of rotation about the axis which is. the, 
eigenvector of the E. matrix- 

gimbal angle 

value of at which L_^ saturates 
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W’a 

W 9 b 

l *'c’ , *’c’ e c 


CO . 

ai 


Euler angles -of the actual attitude of MADD 

Euler angles of the MSS body fixed axes 

Euler angles for the command attitude coordinate 
f rame. 

actual MADD body angular jrates relative to the 
moving; coordinate frame, (i = x,y,z) 

''MSS: angular rates about principal .axes , (i = l" v 2 ,3) 

gyro angular velocity 
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I. Introduction 

In the operation of any manned space vehicle there is a small but 
finite probability that an accident will render it disabled and tumbling. 
It is probable that such a vehicle must be detumbled before evacuating 
the crew and repairs can be performed. ''Tumbling" is the res.ult of a 
significant attitude perturbation to an uncontrolled* body. This results 
in continuous angular motion about all ■three principal body axes, i.e.» 
no inertially oriented axis. The vehicle would reach a stable spin 
after a sufficient amount of energy is dissipated, but this might take 
many weeks or months for a large craft. Astronauts trapped in such 
a body could not easily escape. 

Situations that are most likely to cause tumbling are collision, 
malfunctioning thruster, and escaping stored gas or liquid. From a 
worst case analysis, estimated tumble rates of a large modular space 
station (MSS) are as high as 2 RPM about the principal axes.'*' Elimina- 
tion of angular motion must be done from the nontumbling frame of the 
shuttle orbiter. A Module for Automatic Dock and Detumble (MADD) would 
perform an orbit transfer from the shuttle and then execute a. docking 
maneuver with the vehicle. Internal devices could be employed to reduce 
tumbling to simple spin. However, MADD is adaptable to various 
situations and can be used for the simplified case (spin) but only the 
general case (tumbling) is treated here. Once docked MADD will apply 
torques by firing its thrusters to de tumble the craft. This will be 
done in a time optimal manner. Then the crew can be rescued and the 
mission completed. 
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II. Description of I-IADD 

The purpose of MADD is to de tumble a large passive space base. The 

MSS will be used as a realistic example with initial tumbling at an 

equivalent rate of about 2.RPM. Conceptual design of MADD was influenced 

by shuttle mission objectives and systems constraints. This module must 

be able to maneuver about, dock with, and detumble the MSS with a limited 

amount of fuel for various tumbling situations. The size of MADD is 

limited by the need for manuverability and shuttle specifications, 

29,500 Kg (65,000 lb) payload and payload' bay 18.3 m (60 ft) long by 

? 

4.5 m (15 ft) in diameter.” The size is also somewhat dependent upon 
the configuration of the target vehicle. 

It is necessary that the. CM of. MADD remains fairly well fixed as 
fu'eT* is" consumed and that* mbmen'ts of" inertia change little during 
transfer to simplify control requirements. Therefore, the fuel tanks 
should be located as to minimize this effect. The body frame should 
coincide with the principal axes to simplify control requirements, and 
thrusters should be far enough from the CM to minimize attitude propellant 
usage during... transfer, but not so far as to put excessive moments on 
the docking mechanism during detumbling. MADD must have full orbit 
and attitude control for transfer to the target. A preliminary configura- 
tion for this device is shown in Figure b. 

The structure of this vehicle contains all subsystems with the 
docking -probe- mounted beneath and docking, drogue mounted above the 
main structure, which is a 2.74 m (9 ft) octagon and is 1.22 m (4 ft) 
deep. This contains propellant (hydrazine) storage tanks, control 
systems, batteries, thrusters, and twin-gyro controllers. Subsystems 
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include structure, digital computer, command and telemetry, power 

supply, control systems, guidance, and docking aids. The power supply 

will consist of storage batteries. Control systems consist of three 

types: position, attitude, and detumble. The position and attitude 

control systems provide maintenance of orientation and position- during 

transfer and docking. Twin-gyro controllers' and thrusters were chosen . 

for attitude control devices. Monopropellant hydrazine thrusters were 

chosen for position control, detumbling, and momentum dumping-, because 

3 

of- handling properties, specific impulse, and proven reliability. 

During docking and detumbling the twin-gyro controllers will be locked. 
Once the disabled vehicle is detumbled, the twin-gyro controllers may be 
given new reference signals and released, and the position control system,, 
may be reactivated. Thrusters used for position control are also used 
for detumbling and attitude control. This was done to eliminate the 
need for separate systems of thrusters, even though there are three 
separate control systems. The thrust profile during the detumbling 
procedure is computed by the on-board computer according to the rules 
of optimum detumbling logic. The guidance system consists of an inertial 
platform that provides information on position and orientation of MADD, 
and rate gyros provide body rate information. 

Docking aids consist of laser radar and corner reflectors and 
these are discussed further in the next section. The docking apparatus 
consists of two separate systems ;■ one for a space base and one for. the 
shuttle. A MADD docking probe is used to dock with one of the docking 
ports of the target vehicle, e.g., one located at the extremities of 
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the MSS modules. Capture latches are activated automatically after 
alignment. The shuttle docking systems consists of a remote-controlled 
manipulator boom, probe, capture latches, and a docking droque, which 
is controlled by shuttle crewmen. MADD is positioned about the shuttle 
as illustrated in Figure 2. It is then released to perform the mission, 
afterward the manipulator boom is used to retrieve MM>D and stow it in 
the payload bay. 
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Rescue operations begin as the shuttle completes its rendezvous 
with the distressed vehicle. The MSS will be used to illustrate the 
procedure's. A typical rendezvous is shown in Figure 3. It is assumed 
that the methods required to locate the MSS, determine angular rates 
and momentum are available to the shuttle crew. A stand-off position 
will be established .approximately 61 m (200 ft) away and along the orbital 
path. Mission profile for MADD after deployment from the cargo- bay can 
be broken down into three phases: (1) thrust-free orbital transfer to a 

rendezvous point; (2) thrusted pacing with the docking port and docking.; 
and (3) detumbling of the MSS. MADD becomes automatic at the rendezvous 
point and data is telemetered to the shuttle. Radio and visual contact 
may be, lost intermittently because of occultation. 

The rendezvous point is chosen so that the velocity vector at the 
instant MADD reaches it, will coincide with the velocity vector in the -MSS 
moving coordinate frame. This would eliminate the need for a terminal 
maneuver to reorient the velocity vector, (they should coincide) and 
als.o. reduce the. risk of a collision with the MSS. Another constraint., is 
that the trajectory must not allow MADD to collide with the MSS on its 
way to the rendezvous point. The actual transfer to the rendezvous 
p6int may require several impulses and corrections. An ideal transfer 
would require only one impulse. At the rendezvous point MADD should 
be approximately 4.6m (15 ft) av T ay from the docking port. Thrusters 
begin firing to keep pace with the docking port maintaining proper 
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position and attitude, while closing into, dock (shown in Figure 4)... 

Automated rendezvous and docking is provided by a laser radar 

located on MADD and corner reflectors located on the docking port of 
4 

the MSS. MADD carries a transceiver that provides range, range rate, 
angles-, and angular rates with respect to the docking port. This- radar 
uses CW. (continuous wave) modulation of : an incoherent- gallium arsenide 
injection diode laser. Acquisition occurs when MADD reaches the rendezvous 
point, at this point the system.will.be turned on and it will recieve 
signals from its own transmitter -reflected -back by corner reflectors 
located on the docking port. This permits proper alignment during 
closure and docking by providing guidance information. MADD continues 
closure until the docking probe has engaged and capture latches are 
secured. 

Qnce docked angular rates are measured, and a thrust profile is, 

computed. Thrusters can then detumble the MSS quickly, and rescue is 
completed. 
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IV.' Transfer Trajectories 

Transfer of MADD from shuttle to MSS is divided into two phases. 

Phase one is idealized as having one applied impulse chosen so that 
the vehicle has appropriate position and velocity components at the 
rendezvous point.. During phase-, two. MADD paces the docking port while 
closing into dock. Thrust computations are performed in a moving 
coordinate frame with the origin at the CM of the MSS as illustrated 
in Figure 3. The X-axis is along the direction of motion, Y-axis 
is' normal to the. orbital plane,, and. Z.-axis is along the local vertical*.- 
Equations of. motion for a transfer trajectory to an object in a nearly 
circular oribt are well-known. ^ In the moving X,Y,Z frame during phase 
one, these equations are 

X •+ 2nz = 0 

y + n 2 y = 0 (1) 

z - 2nx - 3n 2 z = 0 
3 1/2 

where n = (GM^/a ) , the mean motion of the MSS in its orbit. The 

solution of set (1) is readily obtained in closed form: 

2z q 4x 2z 

x(t) = cos nt + ( + 6z ) sin nt + (x - — — ) - (3x + 6nz )t 

n no on o o 

y Q 

y(tj = y Q cos nt + — - sin nt (2) 

• • * 

z 2x 2x 

z(tj = ~ sin nt - ( — — + 3z ) cos nt + ( — — + 4z ) 
n no no 

With initial conditions 

x(0) = x q , y (0) = z (0) = 0 
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-set (2) becomes 


2z Ax 2z q 

x(t) = cos nt + sin nt + (x - ) - 3x t 

n n on o 


/ \ o 

y(t) = — sm nt 
n 

• • 

z 2x 2x 

» * o o , o 

z(t) = — sm nt - cos nt H 

n n n 


( 3 ) 


Initial conditions are based on the assumption that the shuttle 
is in the orbital plane of the MSS. The out-of-plane Y-component 
results in simple harmonic motion, while in-plane transfer motion 

is coupled. The only acceptable values of initial conditions x o> 

* 

x , y , and z are those which result in x, y, z and x, y, z 
o o o 

• • • 

simultaneously approaching the values x^, y^, z^ and x^, y^, z^, 
respectively at time t = t^, (rendezvous point). The initial velocity 
components are given as 


. n(x, - x ) sm nt, + 2nz, (1 - cos nt, ) 
i o i l i 

x = 

8(1 - cos nt^) - 3nt^ sin nt^ 

ny ] _ 

y = — : 

o sm nt-^ 

. nz, (A sin nt, - 3nt,) - 2n(x, - x ) (1 - cos nt,) 
i 1 i 1 o 1 

z = ; 

8(1 - cos nt^) - 3nt^ sin nt^ 


( 4 ) 


These resulting expressions indicate that the initial relative veloci'.ty 
requirements for transfer to the rendezvous point are functions of x q , 
x^, y ^ , z^, and time of transfer, t^. x q is dependent upon safety of 
the shuttle to prevent a hazardous situation, and x^ , y , z^ are dependent 
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t^. The velocity 


(5) 


At time t^ the velocity components of’ MADD should be equal to the 
velocity that point x^, y^, z^ would have if it were fixed to the MSS. 
Therefore, the velocity components at time t^ are^ 

• 

xv =- to z, - to y, = -2z sin. nt,. + 4x cos nt, - 3x 
1 y 1 z.l o 1 o 1 o 

f l- = Vl “ Vl “ y o COg Flt l C6) ' 

• • 

z, = CO v - to x, = z cos nt, + 2x sin nt, 

1 x'-l y 1 c 1 o 1 

From these previous expressions, time t^ and initial velocity requirements 
may be determined. ( 0 ^, ( 0 ^, y^ and z^ are relative to the moving 

coordinate frame. They are related to the body fixed coordinate frame 


Upon the location of the 'docking port at time t = 
components are given as 

x(t) = -2z sin nt + 4x cos nt - 3x 
o o o 

•- , • 

y (t) = y Q cos nt 

• • •• 

z(t) = z cos nt + 2x sin nt 

o o 


of the MSS by the following expressions-. 


X 1 


V 

; ' y l 

Z 1 

= T (t) 

' y bl 


t = t. 


(7) 
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■Mi 

X 


*w 

CO 

y 

= T(t) 

10 

by 

to 

z 


>z; 




where 


- 

C4>, cti - s<j>, s* ce^ 

b b b b b 

- s »„ c \ - c * b % . c% 

se b .,s^ b 


T (t) = 

c<i> b s^ b - sc? b c^ b , ce b 

Sfc s W + C.<j>:_ Cfr C0, 

b ; b' b' b d 

- s6 b c v 

(9) 


se b 

se b c<p h 

ce b 


using C? b = cos 

4> b » S 4> b = sin etc - 





Automatic position control during phase two can be modeled by using 
the nonhomogeneous form of equations (1) . 

x + 2nx = f 

x 

y + n 2 y = f -(10) 

z - 2nx - 3n 2 z = f 

z 

where f^ f , f ^ are the applied acceleration components which are the 
control and disturbance forces. Initial conditions associated with set 
(10) become 

x(0) = x 1? y (0) = y 1 , z(0) = 
x(0) = x.^, y (0) = y^, z (0) = z 1 

Taking the Laplace .transform of the differential equations and solving for 
X(s), Y(s), and Z(s) gives 


X(s) = - x - 


.6a.~ 


2 2 

, s - 3n ’ 2n 

z, + — — r x, - 


s i 2,2 2, '1 2 , 2 2.1 , 2 , 2. ‘l 

s (s + n )■ s (s +- n ) s (s + n ) 
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+ s2 : 3 " 2 . f (s> - — 2n 


, 2 , 2 . x 
s (s + n ) 


, 2,2. z 
s (s + n ) 


F (s) 


Y(s) = 


Sy i + y l 1 

-fr^ + V s) 

s + n s + n 


Z(s) = 


2 2 

1 ' . s + 4n 2n 

z_ + - — z, + 


22 1 / 2 , 2 ‘ 1 / 2 2 * A 1 

s + n s (s + n ) s(s + n ) 


1 F (s) + ■. 0 2n 0 F__(s) 


2 2 z 

s + n 


,22 x 
s (s + n ) 


The control law developed is of the form 


f = - K (x - K x) 
CX X 1 


f = _ K (y - K y) 

cy y 2 


( 12 ) 


f = -K (z -K_z) 
cz z 3 


where x, y, z and x, y, z are position and velocity errors with respect 
to the moving coordinate frame. This control law requires position, 
velocity, and attitude information. It also requires the use of throttle- 
able thrusters or a multifunctional monopropellant propulsion subsystem. 
The negative signs in the brackets assure negative feedback. Control 
values , X c » Y Z are computed by the onboard digital computer--. 

Control values f^, f^* f ^ related to the body fixed coordinate frame 
of MADD by the following expression: ^ 
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“ — 
f 

' 

* - 
f 

cl 
f - 

= T (t) 

cx 

f 

c2 
f . 

cy 

f 

c3 


cz 


where 


( 13 ) 


ULL = 


Cd> Cip - Scb 

c c c 

— Scb Cip - Cd) 

c c c 

S0 Sip 
c c 


ce sij> 

c c 

Cd) Sik •+ -Sd> -G0- Cip - 
c r c c c c 

S<j) 

. c 

C3 Si b 

c c 

-Si Sip + Cd> Cd" Cip" 
Y c ' c c ■ c c 

C(j)- 

c 


- se cty. 

c c 




using C<f> c = cos <j> c> S4> c 


sin d> etc. 
c 


Figures 5, 6, and 7 illustrate the block diagrams for the three 
components relative to the moving coordinate frame. The X and X 
components are coupled and the Y component is uncoupled. System 
transfer functions are determined by taking each component separately. 
The X component has inner loop transfer function. 



2 2 
K (s Z - 3n ) 

-■--S 

, 3. . „ 2 , 2 . v .2* 

s(s + Ks + n s - 3K n ) 
x x 


GH = 


2 2 

K (s Z - 3n ) 
x 

•2 2 

s (s + n ) 


(15) 


and outer loop transfer function 


K v K 1 (s 2 - 3n 2 ) 


3 2 2 2 

J ' "i K, + r“'“ - 17 - - 

X X 1 X 


X c s^+Ks J + (K K + n~)s^ - sK n^s - 3^ K n 

' x 1 


( 16 ) 


K .^(s 2 - 3n“) 


3 2 2 9 

s (s + K s +ns-3K n“) 


GH = 
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Figure 5. Block Diagram of X and Z Coupled Position Control System 
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The y component has inner loop, trans.fer function 




€2 2 
2 s z + k s + n 
X 


K s " 

GH = 


2,2 
s +, n 


and outer loop transfer function 


Y _ K v K 2 

Y s 2 + K s + (K K. + n 2 ) 

c y y 2 


GH = 


K K„ 

v2 


2 . „ ,2 

s -r k s -r n 

y 


The Z component has inner loop transfer function 


K 


£ 3 s 2 + K s + n 2 
z 


GH = 


K s 

Z: 


2 , 2 
s + n 


and outer loop transfer function 


Z_ 

Z 


K 

z 3 


s 2 + K s + (K K. + n 2 ) 
2 z 3 


CM = 


K K 
z 3 


2 2 
s + K s + n 
z 


(18) , 


(19) 


( 20 ) 
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Steady state errors of the systems for position, velocity, 
acceleration were -found for each component. For the X control 
(class 1 canonical feedback system) 

. lim s(s 3 + K s 2 + n 2 s - 3K n 2 ) 

, V 1 x x 

e (°°) = n— s-H) x = 0 

P p K K (s 2 - 3n 2 ) 

x 1 


^ lim 

e (°°) - — s-K) 
v K 

v 


3 2 2 " ~ -~2 ~ ~ 

s J + K s + n s - 3K n 
__x x _ 1_ 

K K (s 2 - 3n 2 ) K 1 

X X 


3 2 2 2 

lim s + K s + ns - 3K n 

e (") = — = s-*0 — = 00 

a K K.(s - 3n Z ) S 

x 1 


For the Y control system (class 0 canonical feedback system) 


K = 


lim 

s-K) 


s^ + K s + n Z 

Y. 

K.K 
2 y 


n 


K K 

2 y 


e (oo) = 1 

P 


n 


n 


1 + K 2 K K Vi 'y '2 

p n + K y K 2 y 2 * 

^ lim s 2 + K s + n 2 

p (co) = — = S"K3 — — = oo 

v- K K K 

v y 2 

^ lim s 2 + K s + n 2 
e (°°) = rr~ = s-K) ^ = °° 

3 . K. T - Z 

a K K s 

y 2 


(K K. » n ) 


and 

system 
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For the Z control 'system (class 0 canonical feedback system) 

2 o 

S (“) = — 2 — (K K >> n Z ) 

P’ K K 2 3s 

z 3 


g (oo) =; CO 

V 


g (oo) = oo 

.• a 


The gain of the coupling terms for the X and Z components are 
of the order of the mean motion n, the coupling has relatively l-ittl-e 
effect on the x and z systems. Therefore, the coupling term will be 
neglected temporarily. 

The X component function given by equation (11) is complicated 

by the . two zeros in the numerator of the system equation, (+ /3n. 

Considering the magnitude of . the mean motion, an apporximation is 

2 2 2 

made by cancelling the s - 3n term with s in the demoninator. 

These simplications result in three uncoupled second-order component 


systems . 

The root locus technique was used to determine the performance 
of the component system. ^ Coupling effects and initial conditions 
w.ere neglected for the stability analysis. Since coupling is so 
small, stability should be effected little. The inner and outer 
loop root locus- diagrams ■ for the nOhslmpli-fied X- component systems' 
are shown in Figures 8 and- 9, respectively. The inner and outer 
root locus diagrams for the y component -system are shown in Figures 
10 and 11, respectively. The Z component root locus are similar to 


the Y root locus. 
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Figure 8. Root Locus Plot for X' Component (inner loop) 
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Figure 9. Root: Locus Plot for X Componervt (outer loop) 
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.Since .approximations involving, the mean mot.io.n..and coupling effects 

between the X and Z components resulted in similar transfer functions, 

the, stability analysis for each system will be similar to that- of' the 

2 2 2 

Y system. Validity of cancelling the s - 3n term with the s on 
stability of the X component system can be shown qualitatively. 

The original transfer function has three poles and two zeros in the 
vicinity of the origin, poles at ^ K^/2, + /3n, and zero at + ^3n . 

One zero and one pole .(+ /3n) lie in the right-half plane; indicating 
instability of the system. From the Routh's table the system was- 
found to be unstable for any value of gain. However, the instability 
is small and will respond long before an instability can reach a significant 
value . 

The root locus for an example control system for the Y - component 

is shown in Figure 12. The gains and were chosen to be 100 and 80 

-3 

respectively and n is given as 1.1 x 10 rad/sec. The characteristic 
equation for a second order system is given by 

s 2 + 2£w s + w 2 = s 2 + K s + (K K + n 2 ) (21) 

n. n y y L 

Substituting in values of and K^, the conditional frequency, co = 

80 rad/sec, the natural frequency, to = 89.5 rad/sec, and the damping 
ratio, £ =0.45. The maximum overshoot for a unit step function 
input is given- by. 


- tt £/l - £ -1.57 

e = e 


0.20 


M = 1 + 0.20 = 1.20 
P 
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cus 
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The rise time is given by the equation 


T = 
’ r 


“lO t - "90 t 

CO 


_ = 2,0 - .0.4 _ 1.60 

r 89.5 89.5 


T =~0.018 sec 

r v .. - 


A digital simulation of system response was made in order to check 
the validity of approximations made in the analysis. The digital 
simulation was simplified by having f^ = 0 and initial conditions,, 
zero. The resultant block diagram is shown in Figure 13. Gain 
constants used were the ones used in the previous analysis, = 80 
r and = 100. The Y component system response in position was found 
for unit step, ramp, and sine inputs. The sine input had a period 
of 30 sec. (2RPM) and approximates Y c for a tumbling situation. The 
response for these inputs are shown in Figure 14. Significance of the 
mean motion mangitude on response was found to be negligible for these 
values of gain constants. For a sine input (f = 2 RPM) with magnitude 
of 18.3 m (60 ft) the error was approximately- 3 cm (1.2 in). Thus, 
the response of the Y-component system represents all three systems. 

From these considerations it seems possible to use three simple 
control systems as illustrated in Figure 7 to maintain the proper 
position of MADD during phase two of the mission. 
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Figure 14. Responses for Simplified Y Component Position Control System 
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V. Attitude Control Systems 

Development of the attitude control systems for MADD are considered 
in- this section. Two attitude control systems are required during the 
mission: (1) momentum exchange during phase one; and (2) mass explusion 

during phase two. Attitude control is not required during phase three 
since MADD is docked with the MSS and being detuiiibled. MADD stabilization 
and control systems (SCS) normal modes of operation, include: (1) 

Acquisition establishes the desired three-axis orientation upon deploy- 
ment of MADD from the shuttle; (2) Reorientation includes the slew 
and capture of MADD from one known orientation to another with continuous, 
attitude during the maneuvers; and (3) AV mode provides the SCS configura- 
tion with continuous attitude control during periods of velocity change. 

Requirements and constraints on the control systems for MADD include 
the following areas: (1) sensor type, thruster location, disturbances, 

MADD configuration; (2) stabilization during velocity - change (AV) 
maneuvers; (3) stabilization and control accuracy and their leverage on 
system design; (4) attitude rate control; (5) attitude maneuverability; 

(6) three control degrees of freedom; (7) SCS modes of control; (8) 

MADD dynamic model; (9) MADD static model; and (10) weight restrictions, 
power availability, and reliability. 

Utilization of software techniques requires accurate dynamic 
models in the following areas: (1) accurate rigid-body dynamics, 

including center-of-mass shifts, inertia values, and physical locations; 
(2) sensor models, including dynamic and stochastic error representations 
of sensor, noise; and (3) torque-producing me chanisim dynamics ,. including 
nonlinearities and realistic nonrepeatabiltiy models (stochastic parameter 
variations). This model accuracy is needed for the application of 
sophisticated onboard data-p rocessing techniques. Because of the very 
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short duration of the mission elastic-body dynamics and nonstatic 

physical modeling effects such as thermal deformation are small and 

will not be needed in the dynamic model of MADD. 

Twin gyro controllers and thrusters are used for attitude control. 

Twin gyro controllers are used for phase one and have the following 

advantages: (1) first order .cross-coupling. terms are -eliminated by 

using two-counter -rotating gyros; (2) less, power., and. weight are required 

for a given momentum exchange capability; and (3) larger gimbal angles 

may be used so that a major portion of the stored momentum can be 

transferred to MADD. The equations of motion for a twin gyro controller 
9 

are 

I p = - 2 C fl 6 + L 

X r Z Z Z X 

I q = - 2 C 2 5 4- L. (22) 

y n x x x y 

I r = - 2 C ft 6 + L 

z y y y z 

The small perturbation approach was used to uncouple the equations, thus, 
the second order terms can be neglected. From the equations, it is 
apparent that the controller on the X-axis controls the Y-axis-,. Y controls 
Z, and Z controls X, respectively. By using identical gyros on all three 
axes . 

C = C = C = C and ft = ft = ft = ft 
x y z x- y z 

Therefore, 1 the control equations for each axis are the same except for 
the moment of inertia about that axis. Hie three axes can be- controlled 
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by *an identical control system except for the system gain whicji would 

vary with the moments of inertia. A block diagram for the yaw rate 

g 

system is shown in Figure 15 . The root-locus of this control system 
is shown in Figure 16. 

Attitude control for phase two is accomplished by using thrusters 

placed about the control axes. Mass-expulsion devices are . generally 

inefficient as compared to momentum-exchange devices. Thrusters were 

chosen because the spin vector in not inertially oriented, and momentum- 

exchange devices are incapable of continuously reorienting the spin 

vector without continuous momentum dumping. The attitude control system 

must allow the simultaneous rotation about all three body axes of MADD. 

The control law is symthesized from Euler's theorem on rotation which 

states that any attitude change of a rigid body, may be accomplished- 

by a single rotation about a properly chosen axis. This axis of rotation 

is the eigenvector of the direction cosine matrix defining the orientation 

-of MADD's present position with respect to the command position. 

A control law must be specified that will drive the present attitude 

•of MADD, given by a 3 x 3 direction cosine matrix defined as A to the 

— £ S 

command attitude., given by a 3 x 3 direction cosine matrix defined as 

A in a continuous manner. 

- cs 

Orientation of MADD and command orientation are expressed’ in terms 
of the moving coordinate system by 
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Figure 15. Yaw Rate Control System 
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-1 T' 

-Since A is an -orthogonal matrix (A = A ) and from equation 14 

C S- C S _ CS-v . ' 

the- orientation of MADD with respect to the command orientation is 
gi-ven by 


— 

, 

■_ - 

a. 


c. 

_1 

T 

_1 

a„ 

= A A 

C„ 

_2 

. as cs 

2 

a- 


c_ 

•3' 


3 

. j 


- m 


(24) 


MADD will attain the command attitude when the matrix product A A 

r as cs 


forms an. identity matrix. The matrix product A A 

as cs 


= E is- 


defined as the error matrix. Control law information will be extracted 
from this matrix. 

The elements- of A will be measured by an onboard inertial 
as 

platform and are given by the following matrix 



. Cep. Cip - Sep. C0 Sip 
a a a a a 

C4) Sip + Sep - C0- 
a a a a 

Sep 

a 

A 

-Sep Cip - Cep C0 Sip 

-Sep Sip + Cep C-.0 

Cep 

as 

a a a a a 

a a a a 

a 


se sip 

a a 

-se Cip 

a r a 

ee 

a 


(25) 


using Cep. = cos <p , Sip = sin <p , etc. 

3. a a a 

Where ip^, <p^.> and 0^ are the measured Euler angles of MADD with respect 
to the moving frame. The geometry of the coordinate frames" relative 
to the-rcoving frame are shown in Figure 17. The coordinate frame of 
MADD' is translated in the Z direction and the X' and'- Y axes remain 
parallel to those of the docking reference. The docking reference 
frame is translated relative to the body fixed axes of the MSS and 
is fixed. With these conditions, proper orientation of MADD will 
occur when Euler angles and rates are equal to those of the body fixed 
frame relative to the moving coordinate frame. Since all the elements 
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•of -A are known and the elements of A can be cpmputed because of the 
previous restrictions, the error matrix E_ can be 'computed from Equation 
2 . 4 . 

According to Euler's theorem on rotation, the eigenvectors of the 
matrix b, is the axis about which MADD is rotated to obtain the 
command orientation. The angle of rotation about this axis- is 3. The 
control law is developed such that control torques to the torques L^, 

L , and L do not exceed design values but the system response must 

y z 


be fast. The control law developed is of the form' 
1 . L (max) 


10 


L. = - 

l 


w . 

ai 


I (max) 


[SAT (3, 3 ) b. + , v 

s i 0) (max) 


] 


( 26 ) 


(i = x,y,z) 


L(max), I (max), co o (max) are the maximum values of control torques, 

moments of inertia, and body rate respectively. Ther term SAT (3,. 3 .) 

s 

is a given term derived from the error matrix E: 


SAT (3, 3 ) = 1 for 3 > 3 
s — s 

= 3/3 for 3 < 3 C 

b S 


( 27 ) 


The system saturation level is set by this term. Since L(max), I (max), 
oo (max) are constants for a given rescue mission the control law can 
be simplified' to give 


L 


i 



(i = x,y,z) 


( 28 ) 
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.The terms of b^, b^, and b^ are the components of the eigenvector b in 
the respective body axes and are given by 


1 f E 23 ~ E 32' = 1 E 31 E 13 

•2 -sin 6 ) y 2 sin 6 


b 

z 


1 E 12 ~ E 21 
2 sin 8 


) 


where 


E = 


E 11 

E 12 

E 13 

E 21 

E 22 

E 23 

E 3* 

E 32 

E 33 


(29) 


(30) 


Because of disturbance torques and nonsymmetrical mass configuration 
the components of the eigenvector are time-variant. 

The angle 8 is derived from the. E_ matrix and is given by 


8 = cos 1 [|(E n + E 22 + E 33 - 1) ] (31) 

A block diagram of the attitude control system is shown in Figure 
18. The attitude information from the inertial platform and rate gyros 
is converted from analog form to digital form by the digital computer 
by an analog to digital converter (ADC). Digital orientation informa- 
tion is processed to form the A matrix. This is multiplied with the 

ciS 

T 

command A^.. matrix . and the -error matrix E is generated. - The.- components 
of' the" eigenvector b' and' 8 are computed and body' rates are -suppli'eAT 
to generate the control law. Control law information is converted 
from digital to analog, form by a digital to analog converter (DAC) . 
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Figure 13. Block Diagram of Attitude Control System 
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This will provide, the thrusters with, a continuous thrust profile to. 
orient MADD during phase two of the mission. 

A digital simulation of the attitude control system was made to 
check the accuracy of the control law, (Equation 26), for a given 
tumbling situation. The command attitude was computed by using the 
rigid body euler moment equations with zero' moment applied -to ‘the MSS;. 
The actual attitude and angular rates of MADD were computed by using 
the rigid body euler moment equations with control moments applied. to 
MAD'D". 

For the simulation the following assumptions were made and initial 
conditions chosen. For the- MSS at t = 0 


i* 1 = 0.01 RPM, uj 2 = 0.00 RPM, U 3 = 1.00 RPM 


\p = 0.0°, <p = 0.0°, 0 = 0 T 


where 


0 L = sin 1 


- * 

T - I 


13 


1 - I 


13 


1/2 


I =. — * + — — p - — 

13 A* T ’ max : 2G 

max 


T- = ~(A co i 2 + B <v 2 2 + C oo 3 2 ) 


H 2 = A 2 ,u 1 2 + B 2 w 2 2 + C 2 « 3 2 


For MADD at t =0 

Uj. = 0.01 REM, U)^ = 0.00 RPM, to 3 = 1.00 RPM 

Ip = 0 . 10 °, (p =- 0 . 10 °, 0 = 0 + 0 . 10 ° 

L 
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II = 1 2 = 5.6 .2 Kg-m 2 (1..33 s x 10 3 lbjn-ft^, I 3 = 423 Kg-m 2 a. 00 x lp. 4 Ibm-f.t 2 ) 

cj (max) = 1.00 RPM 
a 

6g = 10° 

L(max) = (0. X ft-lb) _ _ _ - 

The ; rendezvous point was chosen to be 

X = - 2 .17m(-7. 13 ft),. Y = - .915m(-3.00 ft), Z = 24.2m(79.5 f-t) 

(MSS body fixed frame) 

closing rate during docking was 3.04cm/sec(0. 1 ft/sec). 

The responses of MADD's attitude control system are shown in Figure 
19. Part (a) shows that angular rates and part (b) shows the euler 
angles during phase two of the mission. The actual euler angles had a 
maximum error of 0.5° relative to the' command attitude during closure 
and the accuracy of the euler angles were found to be controlled by 
the proper choice of L(max). From this simulation it seems that the 
control law and technique for applying this control law is capable 
of maintaining the proper attitude of MADD during phase two of the 
mission. 

The- control-torquer .assemblies during phase two consist of: 

1) propulsion and 2) motors and actuators, etc. MADD attitude-control 
system requires the capability of a single-propellant supply system 
to provide precise continuous thrust levels within a given thrust range. 
Different attitude-control thrust levels are achieved by either a throttleable 




1 System 
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thruster or a multifunctional monopropellant propulsion subsystem. 

.The propulsion system will use monopropellant hydrazine for position 
and attitude control. The throttleable thruster will have one-to-two- 
orders-of-magnitude thrust variability and the multifunctional mono- 
.propellant propulsion subsystem has a f our-to-f ive-orders-of-magnitude 
thrust variability. The choice between the two schemes will depend 
upon thrust requirements as well as desired flexibility between 
different rescue situations. 

Backup modes may use different sensors and/or different control 
schemes than those used to change attitude during operation in the 
normal mode. The backup modes will exhibit some objectionable 
characteristics such as reduced performance, degraded reliability, 
or increased energy expenditure. 

The point may be brought up that since the maximum duration of 
the mission is only a few hours, the need for backup SCS modes may be 
.unnecessary (unwarranted). The reliability of the primary (normal) 
modes might be high enough to make the backup modes redundant and 
penalizing MADD with additional weight and space, unnecessary complexity, 
and increased cost. Tradeoffs maybe made between increased cost of 

primary with increased reliability without a backup and decreased. cost 

■ / 

of primary with backup adding to the total cost. The final choice will 
have to be made depending on the state-of-the-art at the time" of the 
final design of MADD. 

The costs of SCS for MADD will be less with the incorporation of 
general-purpose digital computation* i.e., system standardization. 
Development costs may not be reduced because it is essential that 
system development stay current. Development cost for MADD may be 
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.VI. .Optimal Detumbling Sequences 

The minimum time optimal; detumbling analyses of a distressed- space 
vehicle can be divided into the following areas: constraint on the 

magnitude of control moment vector and constraint on the magnitude of 
each, component of this vector. These areas may be expressed as 


/ 2 2 , 2 . 1 / 2 - 

,(U + u + u 3 ) <_ M 


K 1 - \ 


k = 1,2,3. 


The space vehicle is modeled as a rigid body by Euler’s moment equations-: 


T^ = + ^2^3 (C-B) 

T 2 = Bw 2 + u 1 w 3 (A-C) (32) 

T 3 = Co) 3 + 03 1 w 2 (B-A) 


The objective is to reduce the three angular velocity components to 
zero in minimum time. 

The first type of constraint leads to a fairly simple solution. 11 
It turns out that the required orientation of the control moment is 
opposite to the angular momentum vector and its magnitude is the largest 
available from the reaction jets. Writing 

\ = I k\ k = 1,2,3 

where 1^ = moment- of inertia about k axis and placing into Euler's 


moment equations we get 
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• f r) X (t) + u i Ct) 

x (t) = ^2 {t) 3 1 


* / \ v (,t)' U'o^tO 

x-Ct) « « 2 x 3 (t) i 2 


f x-\ x (t) + u.Ct) 
= a 0 x*Ct) x ? ^ u 3 


x 3 (t) = 2 


where 


Vi3- 

a V r 2 ?>- 


VLii 

a 2 = Vl 


h. ' l l 
a 3 = Vl 


The op 


timal control moment is 


x*(t) 


u min time 


(t) * Ux*(t)U 


(34) 


where 


| x* (t) 


* 2 


* 2 


= [x 1 Ct) + ^ 2 


(t) + 


*\ 2 ct)i m 


and 


x *( t ) is the solution of 


x(t) 


= f[x(t);t] - M 


x(t) 

llx(t)|l 


( 33 ) 
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-starting rfronix^(-t o ) =-'?^ , -i = iy2,-3- -to- origin (x^ = 0, i = 1,2,3). 
For our case.- we ge.t,. for x(t): 


x (t) = a 1 x 2 (t) x 3 (t) - 


M x 3 (t) 


[x^U) + x 2 2 (t) + x 3 2 (t)] 1/2 


x 2 (t) = a 2 x 3 (t) x^t) - 


H x 2 (t) 


2 2 2 1 /2 
tx^Ct) + x 2 (t) + x 3 (t)3 7 


(36) 


x 3 (t) = a 3 x (t). x 2 (t) - 


M x 3 (t) 


[ Xl 2 (t)- + x 2 2 (t) + x 3 3 (t)] 1/2 


These equations were applied to the tumbling MSS caused by collision 
with the Mark II Orbiter assuming 100% kinetic energy exchange. ^ The 
principal axis angular velocities at commencement of thrusting were 
chosen- at t- = 120 sec. after collision; here the or 's are fairly 
large so as to give a good test to this optimization technique. The 
angular velocities are 1.150, 1.750 and -0.445 RPM's about 1, 2 anQ 
3 principal axes, respectively. These velocities were brought to near 
zero in' about, 7 minutes with, the application of a maximum control 
moment vector magnitude (M) of 3400 Nm (2500 ft-lbs). Figure 2.0 shows 
a time history of the principal axis angular velocities during applica- 
tion of the optimum control moment. Figure 21 gives a time history of 
the body fixed thrusts. (lbs)> required, at point x = -2.17 m (-7.13 ft}, 
y = -.915 m (-3.0 ft) and z = 21.2 m (69.5 ft) to give the necessary 
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3400 Nm (2500 ft-lb) moment directed opposite to. the angular momentum 

vector. With the application of 10,200 Nm (7500 ft-lbs) , the MSS is 

detumbled in 140 sec. Figures 22 and 23 show the results. The x,. y,. 

and; z coordinates for the thrust application are the same. Assuming 

I = 240 sec, 45.4 Kg (100 lbs) of fuel are used for the 3400 Nm 
sp 

(2500 ft-lbs) case and 56.0 Kg (123 lbs) are used for the 10,200 Nm- . 

(7500 ft-lbs) case. 

The second type of constraint (u^ <_ nr, i = 1,2,3) present more- 
difficulty in determining the optimum minimum time control moment 
sequence. In this case, the analysis is not as easily accomplished; 
the control moment vector is not simply directed opporite to the 
angular momentum- vector . As subsequent analysis.: will show, the magnitudes 
of the components (u^ of the control moment vector (u) will be the 
largest possible (nr) - what will change will be. the direction of- 
thrust (+, -). This change in directions of thrust (switching times) 
is, in fact, the major concern in this type of analysis. The equations 
describing minimum time detumbling for the constraint u^ <_ nr are as 


follows: 


. 11,12 


x Ct) = a 1 x 2 (t) x 3 (t) + u (t) 


x 2 (t) = a 2 x 3 (t) x 1 ^ t ^ + 


(37) 


where 


x 3 (t) = a 3 x ^ (t) x 2 (t) + u 3 (t) 


x (t) =. I m (t) 
k k k 


co^(t) = angular velocity about k^ axis 


I k = moment of inertia about k th axis, 
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Figure 22. 
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Figure 23. Body Fixed Thrust a X= -2.17m, Y= -.915m, Z- 2'1.2m for 10 >; 200 Nm Max. Torque 
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H - 1 + 


p l^) x ( t ) + p ( t ) '( t ) + _ • . 

2 2 ^ c ' + P 3 x 3 ( t ) 


= 1 + 


Pl<t) a i X -’ (t) V £) + P : (t) u l(t , 

■V«v 3 (‘)V)*p 2 ( l ) Vl) .. 

p 3 (t) 03^(0 Mt> + P 3 (tj „ - (a> 


(38) 


which yields 


P k (t) - - |H. 

3x k 


Pl<t> ' '“ 2 X 3 (t) p 2 (t > -« 3 * 2 (t) p ( t , 


Po(t) 


-“l x 3 (t > P^t) - Vi ( t) (£) 


('39) 


p 3 ' l > - -v 2 a> p I<t ) (l) 


and (for min. 


H) 


VO = u 


which gives - 


-in time * ' V*” 


(40) 


U i Ct) ~ V gn {p 2 (t)} 


U 2 (t) " ® 2 sg n (p 9 ( t )} 


U 3 Ct) “ “ m 3 s § n lp,(t)} 


( 41 ) 



As can be seen from the above equations for the control moment components 

u^(t), u^Ct) and u^Ct), the control history will be of a bang-bang type; 

this requires investigation for the switching times. First order 

13 

gradient techniques may be used to obtain the switching times. 

However, solutions' were not obtained since variable thrusters are 
available on MADD and the mehtod discussed previously. (constraint- . . 
on .the magnitude of , the moment vector) should require less fuel to 
perform an equivalent task. It should be noted that the control 
history will still be bang-bang even if nonprincipal axis are used 
since the equations of motion will remain linear in u^. 
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VII. Conclusions and Recommendations 

Problems related to docking with and detumbling a passive modular 
space station have been considered here. A MADD concept is proposed as 
a means to apply controlled external torques to detumble a manned space 
base without endangering the shuttle and a preliminary design is 
presented. An operational procedure has been outlined and each subsystem 
discussed. Appropriate assumptions on mission- requirements and constraints 
were formulated based on expected future programs and developments. 
Position, attitude, and detumble control systems were developed. As an 
application the MSS is shown to be detumbled in a few minutes, with very 
small reaction jets. Furthermore, it can be stated that structural 
limitations of the MSS and Human tolerance are not exceeded by the 
induced "g" loads, and that the fuel weight is low. The operation of MADD 
has automatic functions, but the shuttle has command control. The. general, 
case (tumbling) has been treated here, but MADD could be used for the 


spin case as well. 
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